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Booth and Cox proposed the E(s 2 ) criterion for constructing two- 
level supersaturated designs. Nguyen [Technometrics 38 (1996) 69- 
73] and Tang and Wu [Canad. J. Statist 25 (1997) 191-201] inde- 
pendently derived a lower bound for E(s 2 ). This lower bound can be 
achieved only when m is a multiple of N — 1, where m is the num- 
ber of factors and is the run size. We present a method that uses 
difference families to construct designs that satisfy this lower bound. 
We also derive better lower bounds for the case where the Nguyen- 
Tang-Wu bound is not achievable. Our bounds cover more cases than 
a bound recently obtained by Butler, Mead, Eskridge and Gilmour 
[J. R. Stat. Soc. Ser. B Stat. Methodol. 63 (2001) 621-632]. New 
_E(s 2 )-optimal designs are obtained by using a computer to search 
for designs that achieve the improved bounds. 

1. Introduction. In an experiment involving m two-level factors, at least 
m + 1 runs are required to estimate all the main effects. A design is called 
supersaturated if the run size is less than m + 1. Under the assumption of 
effect sparsity that only a small number of factors are active, a supersat- 
urated design can provide considerable cost saving in factor screening. Re- 
cently there have been quite a few articles on the analysis and construction of 
such designs. In particular, the E(s 2 ) criterion proposed by Booth and Cox 
(1962) for constructing two-level supersaturated designs was studied by, for 
example, Lin (1993, 1995), Wu (1993), Nguyen (1996), Tang and Wu (1997), 
Cheng (1997), Li and Wu (1997), Butler, Mead, Eskridge and Gilmour (2001), 
Eskridge, Gilmour, Mead, Butler and Travnicek (2001) and Liu and Dean 
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(2002). This article contains further results on the construction of E(s 2 )- 
optimal designs. 

We represent an iV-run supersaturated design for m two-level factors by 
an N x m matrix X of l's and — l's. Each column corresponds to one factor 
and each row defines a factor-level combination. It is essential that no two 
columns of X are completely aliased, that is, there are no two columns 
x and y such that x = y or x = — y. Throughout this article, we also assume 
that each column of X contains the same number of l's and —l's. Therefore, 
N is even and the number of possible factors that can be accommodated is 
at most | ( ) = ( -\ ) • Thus we have 

(1.1) iY-Km<(^-_ 1 1 

The E(s 2 ) criterion, in seeking a design as close to orthogonal as possible, 
minimizes E(s 2 ) = J2i<j s ij I i™) > wnere s ij is t ne (hj)th. entry of X T X. 
Nguyen (1996) and Tang and Wu (1997) independently derived the following 
lower bound for any supersaturated design with m factors and N runs: 

„. , v m — N + 1 o 

(1.2 E(s 2 ) > - — -N . 

V ' V ' ~ (m- l)(N - 1) 

When N = (mod 4) , this bound can be achieved only if m is a multiple 
of N — 1; when N = 2 (mod 4), m needs to be an even multiple of N — 1. 
One question is whether this bound can be achieved for every multiple of 
N — 1 when N = (mod 4) and every even multiple of iV — 1 when N = 2 
(mod 4). This appears to be very hard and the answer is yet unknown. One 
objective of this article is to provide some results in this direction. 

In Section 2 we present a method for constructing designs that achieve the 
Nguyen-Tang- Wu bound. Section 3 contains the other main result of this ar- 
ticle: improved lower bounds for E(s 2 ) when the Nguyen-Tang- Wu bound is 
not achievable. It came to our attention that Butler, Mead, Eskridge and Gilmour 
(2001) also derived improved lower bounds for E(s 2 ). Unlike our bounds, 
Butler, Mead, Eskridge and Gilmour's bounds do not apply to all cases; 
see the discussion in Section 3. We also report some new £'(s 2 )-optimal de- 
signs obtained by using a computer to search for designs that achieve the 
improved bounds. All the proofs are presented in Section 4. 

Throughout this article, we use GF(s) to denote a finite field with s 
elements. The multiplicative group that consists of the nonzero elements of 
GF(s) is cyclic, and a generator of this group is called a primitive element 
of the field. For each positive integer q, we denote the set {0, 1, . . . , q — 1} 
of nonnegative integers less than q by Z 9 . Multiplication and addition in Z g 
are reduced modulo q when necessary. For each subset T = {a\, . . . ,at} of 
Z q and b € Z g , the set {a± + b, . . . ,a t + b} is denoted by T + b. 
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2. Construction of £7(s 2 )-optimal supersaturated designs via balanced 
incomplete block designs. Let m = q(N — 1), where q is a positive inte- 
ger. Cheng (1997) showed that a supersaturated design that achieves the 
lower bound in (1.2) is equivalent to a balanced incomplete block design 
with N — 1 treatments and q(N — 1) blocks of size N/2 — 1, abbreviated 
as BIBD(iV — l,q(N — l),N/2 — 1). Without loss of generality, we may as- 
sume that all the entries in the first row of a supersaturated design X 
are equal to 1. Let Z be obtained from X by deleting the first row. Then 
Z = [%'](iv-i)xm can De considered as the treatment-block incidence matrix 
of a binary incomplete block design with N — 1 treatments and m blocks 
of size N/2 — 1, where the ith treatment appears in the jth block if and 
only if Zij = 1. Then X attains the lower bound in (1.2) if and only if Z is 
the treatment-block incidence matrix of a balanced incomplete block design. 
This equivalence of the existence of a BIBD(A r — l,q(N — l),N/2 — 1) and 
that of an N x q(N — 1) X that attains the lower bound in (1.2) extends 
the well-known result that the existence of a BIBD(iV — 1, N — 1, N/2 — 1) is 
equivalent to that of an N x N Hadamard matrix. Note that for X to have 
no completely aliased columns, all the blocks of the corresponding BIBD 
must be distinct. 

The largest BIBD(iV - l,q(N - l),N/2 - 1) with distinct blocks is the 
trivial one consisting of all the ^ ^ 2 _ 1 1 ^ subsets of size N/2 — 1 of the N — 1 
treatments, which corresponds to the supersaturated design with the maxi- 
mum number of factors given in (1.1). So the question raised in the Introduc- 
tion about the existence of designs that achieve the Nguyen-Tang- Wu bound 
is equivalent to whether there exists a BIBD(iV — 1, m, N/2 — 1) with distinct 
blocks for every m satisfying (1.1) that is a multiple of N — 1 when N = 
(mod 4) and an even multiple of N — 1 when N = 2 (mod 4) . If there exists a 
Hadamard matrix of order N, then a BIBD(A r - 1, q(N - 1), N/2 - 1) exists 
for every q; a simple construction is to piece together q BIBD(iV — 1, N — 
1, N/2 — l)'s. The additional challenge in the supersaturated design construc- 
tion is not to use the same block more than once. In fact, the construction 
of supersaturated designs by combining q Hadamard matrices as proposed 
in Tang and Wu (1997) is equivalent to piecing together q BIBD(iV — 1,N — 
l,N/2 — l)'s. It is not clear how one can avoid duplicated blocks in their 
construction. Also, Tang and Wu's construction is applicable only when 
is a multiple of 4. Furthermore, what one needs is that the whole design is a 
BIBD; it does not have to be the union of q BIBD(iV - 1, N- 1, N/2 - l)'s. 
Thus the construction based on BIBDs is more general and flexible. In this 
section we present a method of using difference families [Wilson (1972)] to 
construct BIBD(iV - l,q(N - l),N/2 - l)'s with distinct blocks. The read- 
ers are referred to Chapter VII of Beth, Jungnickel and Lenz (1999) for a 
discussion of the construction of BIBDs based on difference families. 

The following is our first construction result. 
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Theorem 2.1. Suppose N — 1 is an odd prime power, q is an even 
divisor of N — 2, x is a primitive element of GF(A r — 1) and T is a subset of 
7L q ofsizeq/2. Then the q(N - 1) sets {S r>a :r = 0, . . . ,q- 1, a G GF(2V- 1)}, 
where 5 r>a = {x jq+l + a : < j < (N — 2)/q — 1, i6T + r}, /orm i/ie blocks 
of aBIBb(N -l,q(N -1), N/2-1). Furthermore, if(N-2)/q is odd and 
U is a subset of 7L q of size q/2 such that U* = U + (q/2), where U* is the 
complement of U in Z 9 , then the q(N — l)/2 sets {SV,a - r EU, a G GF(iV — 
1)} form the blocks of a BIBD(iV — l,q(N - 1) /2, N/2 - 1). 

Note that in the above construction, as defined earlier, the elements of 
T + r are reduced modulo q if necessary. 

The BIBDs constructed in Theorem 2.1 may not have distinct blocks. For 
example, if T + r = T, then S ra = <So,a- To construct BIBDs with distinct 
blocks, let e be the smallest positive integer such that T + e = T '. Then e, 
called the order of T, is a divisor of q. Let Tj = {T + (j — l)e, ... ,T + (J — 
l)e + e — 1} for 1 < j < q/e. Then T\ = ■ ■ ■ = T q j e and all the e sets in each 
Tj are distinct. In this case, it can be seen that the design constructed in 
Theorem 2.1 consists of q/e replications of an identical BIBD with distinct 
blocks as long as q ^ N — 2. This leads to the following theorem. 

Theorem 2.2. Suppose N — 1 is an odd prime power, q is an even 
divisor of N — 2 with q ^ N — 2, x is a primitive element of GF(N — 1) and 
T is a subset ofZ q of size q/2. Let e be the smallest positive integer such that 
T + e = T. Then the e(N - 1) sets {S r , a : r = 0, . . . , e - 1, a G GF(iV - 1)}, 
where SV,a = {x^ q+l + a : < j < (N — 2)/q — 1, i G T + r}, are distinct and 
constitute a BIBD(iV - l,e(N - 1), N/2 - 1). Furthermore, if (N - 2)/q is 
odd and U is a subset of size e/2 of {0, . . . , e — 1} such that U* = U + (q/2), 
where U* is the complement ofU in {0, . . . , e — 1} and the addition is reduced 
modulo q, then the e(N — l)/2 sets {<SV,a - r G U, a G GF(iV — 1)} constitute 
a BIBD with distinct blocks. 

The first design described in Theorem 2.2 is constructed by using e initial 
(base) blocks So^, . . . , 5 e _i 5 o, where SV,o = {£ J " 7+ ' :0 < j < (N — 2)/q — l,i G 
T + r}. The second design uses the e/2 initial blocks S r fi, where r G U. This 
construction is similar to Wilson's construction of balanced incomplete block 
designs as described in Theorem 5.2 of Beth, Jungnickel and Lenz [(1999), 
page 489]. However, Wilson did not consider the constraint of no repeated 
blocks. Also, a divisor q of N — 2 was used to construct BIBDs of block size 
N/2 — 1 in Theorems 2.1 and 2.2, while Wilson used it to construct BIBDs 
of block size (N - 2)/q or (N -2)/q + 1. 

We note that if both T and U consist of all the integers < % < q/2 — 1, 
then e = q and U* = U + (q/2). Thus a BIBD(A r - 1, q(N - l),N/2 - 1) with 
distinct blocks can be constructed for each even divisor q of N — 2 and if, in 
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addition, (N - 2)/q is odd, then a BIBD(iV - l,q(N - l)/2, N/2 - 1) with 
distinct blocks also exists. 

Example 2.1. Let N = 20. Then N — 2 = 18 has two even divisors 
not equal to 18: 2 and 6. Since both 18/2 and 18/6 are odd, Theorem 2.2 
can be used to construct BIBDs with 19, 38, 57 and 114 distinct blocks 
of size 9. The construction is based on the finite field GF(19), which is 
equivalent to Zig. Suppose we choose the primitive element 2. To con- 
struct a BIBD(19, 57, 9), we let the integer q in Theorem 2.2 be 6, T = U = 
{0,1,2} and use the three initial blocks {2°, 2 1 , 2 2 , 2 6 , 2 7 , 2 8 , 2 12 , 2 13 , 2 14 } = 
{1, 2, 4, 7, 14, 9, 11, 3, 6}, {2 1 , 2 2 , 2 3 , 2 7 , 2 8 , 2 9 , 2 13 , 2 14 , 2 15 } = {2, 4, 8, 14, 9, 18, 3, 6, 12} 
and {2 2 ,2 3 , 2 4 , 2 8 , 2 9 , 2 10 , 2 14 , 2 15 , 2 16 } = {4, 8, 16, 9, 18, 17, 6, 12, 5}. Adding 
the integers 0,1,..., 18 (mod 19) to all the elements in the initial blocks 
produces a BIBD with 57 distinct blocks of size 9. Write down the 19 x 57 
treatment-block incidence matrix in which the (i,j)th entry is equal to 1 if 
the ith treatment appears in the jth block and is equal to —1 otherwise. 
Then by adding a row of l's to this treatment-block incidence matrix, one 
obtains an £J(s 2 )-optimal 20-run design for 57 factors. 

Note that once the initial blocks are determined, except for the row of l's, 
the other rows of the corresponding supersaturated design can be developed 
cyclically from an initial row. Eskridge, Gilmour, Mead, Butler and Travnicek 
(2001) and Liu and Dean (2002) also considered cyclic generation of E(s 2 )- 
optimal and nearly optimal supersaturated designs. 

The following result provides more flexibility in the construction of E(s 2 )- 
optimal supersaturated designs. 

Theorem 2.3. Let q be an even divisor of N — 2 such that q ^ N — 2. 
Let T and T' be subsets of size q/2 of 7L q such that T' + a for all ele- 
ments a ofTL q . If d± and c?2 are BIBDs constructed by applying Theorem 2.2 
to T and T' , respectively, then d\ and have no blocks in common; therefore 
their union is also a BIBD with distinct blocks. 

Let T be the set of all subsets of size q/2 of Z g . For any two such subsets 
T and T' , we write T ~ T 1 if there is an element r£Z 9 such that T = 
T' + r. Then clearly "~" is an equivalence relationship. Therefore, J- is 
partitioned into disjoint equivalence classes. One can choose a set T from 
each equivalence class to construct a BIBD according to the method of 
Theorem 2.2. The BIBDs constructed by using T's from different equivalence 
classes have no blocks in common. Therefore, the union of these BIBDs is 
a BIBD with distinct blocks of size N/2 — 1, and can be used to construct 
.E(s 2 )-optimal supersaturated designs that attain bound (1.2). Note that 
the order of each set T is equal to the size of the equivalence class that 
contains T. 



G 



D. A. BULUTOGLU AND C.-S. CHENG 



Example 2.1 (continued). Again consider the case N = 20. More de- 



subsets of size 3 of Zg can be partitioned into four equivalence classes: 



Each set in one of the first three equivalence classes can be used to con- 
struct a BIBD(19,57,9) and a BIBD(19, 114,9) with distinct blocks. By 
using one, two or all three of these equivalence classes, one can construct 
BIBD(19,19t,9)'s with distinct blocks for t = 3,6,9, 12, 15 and 18. The last 
equivalence class can be used to construct a BIBD(19, 19, 9) and a BIBD(19, 38, 9) 
with distinct blocks. Combining these with designs constructed from the first 
three equivalence classes, we obtain BIBD(19, 19i, 9)'s with distinct blocks 
and S(s 2 )-optimal 20-run designs with 19i factors, for 1 < t < 20. 

Note that in the above example, the designs constructed by using the last 
equivalence class are the same as those constructed by choosing q = 2. 

Example 2.2. Let N = 18. Then N — 2 = 16 has three even divisors 
not equal to 16: 2, 4 and 8. Let q = 8. It can be seen that the 70 subsets 
of size 4 of 7*8 can be partitioned into 10 equivalence classes: 8 equivalence 
classes of size 8, 1 equivalence class of size 4 and 1 equivalence class of 
size 2. By applying Theorem 2.2, one can construct BIBD(17, 34i, 8)'s with 
no repeated blocks for t = l (using the equivalence class of size 2), 2 (us- 
ing the equivalence class of size 4) and 4 (using each of the 8 equivalence 
classes of size 8). Combining the BIBDs constructed from different equiva- 
lence classes, one obtains BIBD(17, 34t, 8)'s with distinct blocks for all t's 
such that 1< t < 35. Note that since 18 is not a multiple of 4, the number 
of blocks b of a BIBD(17, b, 8) must be a multiple of 34. 

Example 2.3. For N = 10, applying Theorems 2.2 and 2.3 with q = 4, 
one can construct BIBD(9, 18£, 4)'s with no repeated blocks for t = 1,2,3. 
This is because the six subsets of size 2 of Z4 can be partitioned into two 
equivalence classes of sizes 4 and 2, respectively. In this case, the trivial 



BIBD has (4) = 126 blocks. Taking the complements of BIBD(9, 18i, 4)'s 




with no repeated blocks for t = 1,2 and 3 in the trivial BIBD, we obtain 
BIBD(9, 18£, 4)'s with distinct blocks for i = 4, 5 and 6. This provides a 
complete solution of all BIBD(9, b, 4)'s with no repeated blocks, and thus all 
10-run supersaturated designs that attain the Nguyen-Tang-Wu bound. 



signs can be obtained by using Theorem 




{0,1,2}, {1,2,3}, {2,3,4}, {3,4,5}, {4,5,0}, {5,0,1} 
{0,1,3}, {1,2,4}, {2,3,5}, {3,4,0}, {4,5,1}, {5,0,2} 
{0,1,4}, {1,2,5}, {2,3,0}, {3,4,1}, {4,5,2}, {5,0,3} 
{0,2,4}, {1,3,5}. 
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Sometimes one can also produce designs with distinct blocks by combining 
those constructed by using different even divisors of iV — 2. Let Tj be a subset 
of size qi/2 of Z ?i , i = 1, 2, where c/i and c/2 are even divisors of iV — 2. Then 
by the same argument as in the proof of Theorem 2.2, one can show that 
the designs obtained by applying Theorem 2.2 to T\ and T2 have no blocks 
in common as long as all of the initial blocks of the two designs are different. 

Example 2.4. Let N = 14. Then N — 2 = 12 has three even divisors not 
equal to 12: 2, 4 and 6. As in Example 2.1, by choosing q = 6, one can con- 
struct BIBD(13,26t,6)'s with distinct blocks for £ = 1,3,4,6,7,9,10. Note 
that since 12/6 is even, fewer designs can be constructed here than in Exam- 
ple 2.1. For q = 4, since 12/4 is odd, a BIBD(13, 26, 6) and a BIBD(13, 52, 6) 
with distinct blocks can be constructed by applying Theorem 2.2 to the sub- 
set T = {0, 1} of Z4. It can be seen that these two designs have no blocks 
in common with any of those constructed by using q = 6. It follows that 
one can construct BIBD(13, 26t, 6)'s with distinct blocks for 1< t < 12. As 
in Example 2.1, choosing q = 2 does not produce new designs. 



3. Improved lower bounds for E(s 2 ). 

some improved lower bounds for E(s 2 ). 



The following theorem presents 



Theorem 3.1. Suppose m is a positive integer such that m > N — 1. 
Then there is a unique q such that — 2N + 2 < m — q(N — 1) < 2N — 2 and 
(m + q)=2 (mod 4) . Let g(q) = (m + q) 2 N - q 2 N 2 - mN 2 . 

1. If N = (mod 4), then 

( g{q) + 2N 2 - AN 



E{s 2 ) > I 



m(m — 1) 
g(q) - 2N 2 + 4N + 4N\m -q(N-l 



when \m - q(N — 1)| < N — 1, 



m(m — 1) 



g(q)+4N 2 -4N 



when N — 1 < \m - q(N - 1)\ <-N -2, 

3 

when \m - q(N - 1)| > -N - 2. 



m(m — 1) 

2. If N = 2 (mod 4) and q is even, then E(s 2 ) > max(/i(g), 4), where 
( g(g) + 2N 2 -4N + 8^ ^ _ _ < ^ _ 
m[m — 1) 

g(q) - 2N 2 + 20N + (4N - 8)\m - q(N — 1)| — 24 



%) = { 



g(q)+4N 2 -4N 
m(m — 1) 



m(m — 1) 

3 

when N - 1 < \m - q{N - 1)| < -N - 3, 

3 

when \m - q(N - 1)| > -N - 3. 
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3. IfN 



2 (mod 4) and q is odd, then E(s 2 ) > max(/i(g), 4), where 





h(q) 



g(q) - 2N 2 + AN + AN\m - q(N - 1)| 
to(to — 1) 

3 

when N - 1 < \m - q(N - 1)| < -N - 1 

g{q) + AN 2 - 12N + 8\m - q(N - 1) | + 8 
m{m — 1) 

3 

when \m - q(N - 1)| > -N - 1. 



Butler, Mead, Eskridge and Gilmour (2001) also derived some lower bounds 
for E(s 2 ). Write to as m = q'(N — 1) + r, where |r| < N/2. Their result does 
not apply to the case where N = 2 (mod 4) and q' is odd, while our bounds 
apply to all cases. A numerical comparison suggests that their bounds agree 
with ours in the cases where they are applicable. Table 1 shows values 
of the Nguyen-Tang-Wu bound and our improved bound for the range 
N < to < 2(N - 1), where N = 10, 12, 14 and 16. 

We have used the computer software Gendex developed by Nguyen (1996) 
to search for £'(s 2 )-optimal designs. In many cases, we were able to find 
designs which achieve the improved bounds. Since 

citetr3 have reported £?(s 2 )-optimal 12- and 16-run designs, we list in Table 2 
the new 10- and 14-run £'(s 2 )-optimal designs we have found. 

4. Proofs. 



Table 1 

Nguyen-Tang-Wu bound and the bound of Theorem 3.1 
for N <m< 2(N - 1), N = 10, 12, 14 and 16 



N 



rn 



Bound of 
Theorem 3.1 



Nguyen-Tang-Wu 
bound 



10 



10 

11 

12 
13 
14 
15 
16 
17 
18 



4 
4 
4 



1.23456 
2.22222 
3.03030 
3.70370 
4.27350 
4.76190 
5.18518 
5.55555 
5.88235 



4.61538 
5.05494 
5.52381 
5.86666 
5.88235 
5.88235 
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Table 1 
Continued 







Bound of 


Nguyen-Tang- Wu 


N 


m 


Theorem 3.1 


bound 


1 9 
1 Z 


1 9 


9 18181 
Z.IOIOI 


1 i anns 
i.iyuuo 




1 3 




9 18181 

Z. lolOl 




1 1 

1-1 


4 91 Q78 


Q n9flQ7 
o.uzuy i 




1 % 

10 


/I P.7149 


o. / 4UZO 




1 R 


c; 9nnnn 

O.iUUUU 






1 7 


U.U4 l uo 


4 QDOflQ 




1 8 
lo 


o.you I o 


o.oyuo i 




1 Q 

i y 




c. 81 81 8 
o .0 10 10 




on 
_ ■ u 


« S9i nc. 


u.zuyu f 




9 1 
Z 1 


8^.71 /I 


0.04040 






fi 8^71 4 

U.OO / 14 


fi 8^71 4 

U.OU / 14 


14 


14 


4 


1.15976 




i r 

15 


4 


2.15384 




1 (\ 

1U 


A 

-± 


Q ni ^8 
o.uiooo 




1 t 


A 0-41 1 7 


Q 7RQ9Q 

o. / uyzo 




1 X 

lo 


O.U ( ozu 


4 4^4^8 




1 Q 

i y 




c. r)9 c ifi4 




on 
zu 


ft QCi7QQ 

u.ou f oy 






9 1 

Z 1 


ft ftftftftft 
u.uuuuu 


P, D' : in7R 
u.uou / u 




99 
zz 


ft QOOOQ 


fi 4R1 "i^ 




9*3 
zo 


7 41 ^09 


fi S^l 4 




9 1 

Z-l 


7 &9fi0& 
l .ozuuo 


7 9 1 fl7fl 
J .z iu / u 




9^ 

Z-.) 


7 84000 


7 c i'}84fi 




26 


7.84000 


7.84000 




27 


8.38746 


7.87692 




28 


8.80423 


8.21728 




29 


8.82758 


8.53333 




30 


8.82758 


8.82758 


16 


16 


2.13333 


1.13777 




17 


3.76470 


2.13333 




18 


4.18300 


3.01176 




19 


4.49122 


3.79259 




20 


5.38947 


4.49122 




21 


6.09523 


5.12000 




22 


6.64935 


5.68888 




23 


7.08300 


6.20606 




24 


7.42029 


6.78261 




25 


7.68000 


7.11111 




26 


7.87692 


7.50933 



Proof of Theorem 2.1. Since x is a primitive element of GF(iV — 1), 
we have x N ~ 2 = 1, x^ N ~ 2 ^ 2 = — 1 and 1, x, x 2 , . . . , x N ~ 3 are all distinct, 
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Table 2 

New 10- and lA-run E(s 2 )-optimal designs 







N 




10, 


m = 14, E(s 2 ) = 5.0549 












-1 


1 


1 


1 


-1 


1 -1 -1 


-1 -1 1 


-1 


1 


-1 








1 


1 


-1 


1 


1 


-1 -1 -1 


1 1 1 


-1 - 


1 


1 








-1 


-1 


-1 


1 


1 


-1 -1 1 


-1 -1 -1 


1 


1 


1 








-1 


-1 


1 - 


1 


-1 


-1 1 -1 


-1 1 -1 


-1 - 


1 


1 








1 


1 


-1 - 


1 


-1 


-1 1 1 


-1 1 -1 


-1 


1 


-1 








1 


-1 


1 - 


1 


1 


1 -1 1 


-1 -1 1 


-1 - 


1 


-1 








-1 


1 


-1 - 


1 


-1 


1 1 1 


1 -1 1 


1 - 


1 


1 








1 


1 


1 


1 


1 


-1 1 -1 


1 -1 -1 


1 - 


1 


-1 








1 


-1 


1 


1 


-1 


1 1 1 


1 1 1 


1 


1 


1 








-1 


-1 


— 1 - 


1 


1 


1 -1 -1 


1 1 -1 


1 


1 


— 1 














N = 


10, ro = 15, E(s 2 ) = 5.5238 












-1 


-1 


-1 


1 


1 


1 1 1 


-1 1 -1 


-1 - 


1 


-1 


-1 






-1 


-1 


-1 - 


1 


1 


-1 -1 -1 


1 -1 -1 


-1 


1 


1 


1 






1 


1 


-1 


1 


-1 


1 1 -1 


-1 1 -1 


1 


1 


1 


1 






-1 


1 


-1 


1 


-1 


-1 1 1 


1 -1 1 


1 


1 


-1 


1 






1 


1 


1 


1 


1 


-1 -1 1 


-1 -1 -1 


-1 - 


1 


1 


1 






-1 


-1 


1 - 


1 


-1 


1 -1 1 


-1 1 1 


-1 


1 


-1 


1 






1 


-1 


-1 - 


1 


-1 


-1 -1 1 


1 1 1 


1 - 


1 


1 


-1 






1 


1 


1 - 


1 


-1 


1 1 -1 


1 -1 -1 


-1 - 


1 


-1 


-1 






-1 


1 


1 


1 


1 


1 -1 -1 


1 1 1 


1 - 


1 


1 


-1 






1 


-1 


1 - 


1 


1 


-1 1 -1 


-1 -1 1 


1 


1 


-1 


-1 














N - 


= 14, m = 17, E(s 2 ) = 


4.9412 










-1 


1 


1 - 


■1 


-1 


-1 1 1 


1 -1 -1 


1 - 


1 


-1 


-1 


-1 


1 


1 


-1 


-1 


1 


1 


-1 1 -1 


1 -1 1 


1 - 


1 


-1 


1 


1 


-1 


-1 


-1 


1 


1 


1 


1 1 1 


1 -1 -1 


-1 


1 


1 


1 


1 


1 


1 


-1 


-1 - 


1 


1 


-1 1 1 


-1 1 1 


-1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 - 


1 


-1 


1 -1 1 


1 1 -1 


-1 - 


1 


-1 


-1 


1 


-1 


1 


1 


-1 - 


1 


1 


-1 -1 -1 


1 1 -1 


1 


1 


1 


1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 -1 -1 


-1 1 1 


1 - 


1 


1 


-1 


1 


1 


-1 


1 


-1 - 


1 


-1 


1 1 1 


-1 -1 1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


1 


1 


1 -1 1 


1 1 1 


-1 - 


1 


-1 


1 


-1 


1 


-1 


1 


1 - 


1 


1 


-1 -1 -1 


-1 -1 1 


-1 


1 


-1 


-1 


1 


-1 


-1 


-1 


1 - 


1 


1 


1 1 -1 


-1 1 -1 


1 - 


1 


-1 


1 


-1 


-1 


1 


1 


1 


1 


-1 


-1 -1 1 


-1 -1 -1 


-1 - 


1 


1 


1 


-1 


-1 


1 


1 


1 


1 


-1 


1 1 -1 


1 1 1 


-1 


1 


1 


-1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 -1 -1 


-1 -1 -1 


1 


1 


-1 


-1 


-1 


1 
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Table 2 
Continued 











N = 


14, 


m = 


= 18, E(s 2 ) 


= 5.6732 














-1 


1 


-1 


1 


1 


-1 


1 


1 


1 


1 


1 


1 1 - 


■1 


-1 


1 


1 


-1 




1 


1 


-1 


-1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 -1 - 


1 


-1 


1 


-1 


-1 




1 


1 


-1 


-1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 -1 


1 


1 


1 


1 


1 




1 


1 


1 


1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 -1 


1 


-1 


-1 


1 


1 




1 


-1 


1 


1 


-1 


1 


-1 


-1 


1 


1 


1 


1 1 


1 


1 


1 


1 


-1 




-1 


-1 


-1 


1 


-1 


1 


-1 


1 


1 


-1 


-1 


-1 1 - 


■1 


-1 


1 


-1 


1 




-1 


-1 


1 


-1 


1 


-1 


-1 


1 


1 


-1 


1 


-1 -1 - 


■1 


1 


-1 


1 


-1 




-1 


-1 


-1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


1 1 


1 


1 


-1 


-1 


1 




-1 


-1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


-1 


1 1 


1 


-1 


-1 


1 


-1 




-1 


-1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


-1 -1 


1 


1 


1 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 1 - 


■1 


1 


-1 


-1 


1 




1 


-1 


1 


-1 


-1 


-1 


1 


-1 


1 


1 


1 


1 -1 - 


1 


-1 


-1 


-1 


1 




1 


1 


-1 


-1 


_1 


_1 


-1 


1 


-1 


1 


1 


-1 1 


1 


-1 


-1 


-1 


-1 




-1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


-1 -1 - 


1 


1 


1 


-1 


-1 














N-- 


= 14, m 


= 19, E(s 2 


) =6.0585 












-1 


-1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


1 


-1 


1 -1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 


1 


-1 


1 1 - 


1 


1 


-1 


1 


1 


-1 


1 


1 


1 


1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 -1 - 


1 


-1 


-1 


-1 


1 


1 


-1 


-1 


-1 


1 


1 


-1 


1 


-1 


1 


1 


-1 


1 1 - 


1 


-1 


1 


-1 


-1 


-1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


-1 


-1 


1 


1 -1 - 


1 


-1 


1 


-1 


-1 


-1 


-I 


1 


1 


1 


1 


-1 


1 


1 


-1 


-1 


-1 


1 -1 


1 


-1 


-1 


1 


1 


-1 


1 


1 


1 


-1 


-1 


-1 


1 


1 


1 


1 


1 


-1 1 


1 


-1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


-1 


1 


1 


-1 1 - 


1 


1 


-1 


1 


-1 


1 


-1 


1 


1 


-1 


1 


1 


1 


-1 


1 


-1 


1 


-1 1 - 


1 


1 


1 


1 


1 


-1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


1 


-1 


1 


-1 -1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


1 


-1 


1 


1 


1 


-1 


-1 


-1 


-1 -1 - 


1 


1 


1 


1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


1 -1 


1 


1 


-1 


-1 


-1 


-1 


1 


-1 


-1 


1 


1 


1 


-1 


1 


1 


-1 


1 


1 1 


1 


1 


1 


-1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 


-1 


-1 


-1 


-1 


-1 1 


1 


-1 


-1 


-1 


1 


1 



where 1 is the multiplicative identity. The multiset Ur=oLL,ijeS r0 ,M7=«'L u — ^J" 
also can be expressed as 



g-l (N-2)/q-l 

U U U { x j2q+i2+r - x (ji+j2)q+ii+r }, 

H,*aeT r=0 j 2 =0 
0<jl<(N-2)/q-l 
ii7^2 if ji=0 
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where i\ + r and %i + r are reduced modulo q if necessary. For any fixed 
triple (ii,i2,ji), where i\, 12 £ T, < j\ < (N — 2)/q — 1 and *i ^ ^2 if Ji = 0, 

q-l(N-2)/q-l 

U U {x j2,?+i2+r " — x ^ 1+ ^ q+il+r } 

r=0 j 2 =0 

g-1 (N-2)/q-l 

= Q | x i29+*2+r^ _ 2 ,jicj+[(ji+r)-(i2+r)]^|^ 

r=0 j 2 =0 

where again i\+r and 22 + r are reduced modulo (/ if necessary. 

Now Ur=o Uj^Lo 2 ^ 9 1 {x^ 2q+l2+r } covers every power of x and hence every 
nonzero element of GF(iV — 1) exactly once. Since 1 — x^ iq+ ^ n+r '~^ 2+r ^ 7^ 
and there are — 1) triples (ii,i2,ji) such that 21,^2 € T, < ji < 

(N-2)/q-l and n 7^ z 2 if ji = 0, it follows that Ur=o Uu,ves r , ,u^v{ u ~ v } 
covers each nonzero element of GF(iV — 1) f (^T^ ~~ -0 times. Thus the q 
sets in {S r fi : < r < q — 1} are a difference family and the q(N — 1) sets in 
{S r , a :0<r<q-l, a€ GF(iV - 1)} constitute the blocks of a BIBD. 

If (iV — 2)/q is odd, then (q + N — 2)/2 = aq for some positive integer a. 
Then since x^" 2 )/ 2 = -1 and U* = U + (g/2), 

+ r _ x 0'i+i2)<?+ii+»-| 

{ x 9/2 x (A r -2)/2|- x (i 1 +j 2 )g+n+r _ x j 2 q+i 2+r ^ 
X aq ■ 5.x^ 1+ ^ qJril+r — x i2<?+i2+r| 



r ^.(ji+j2+a)q+h+r rr ,(j2+oi)q+i2+r^ 



(TV 


-2)/q-l 


u 




reu 


i2=o 




(JV-2)/ 9 -l 




u u 




ret/* J2=0 




(AT-2)/g-l 




u u 




ret/* J2=0 




(/V_2)/g-l 




u u 




rg{/* jr' 2 =0 




(/V_2)/g-l 




u u 




ret/* J2=0 


This implies 


that Uret/U W) , 



a j(j'l+j2)g+n+r _ x 32q+i2+r\ 



Since Ur=o Uu^e5 r o,«A^ u ~~ v ^ covers eac h nonzero element of GF(iV — 1) 

times, eachof{J reU {J u v£Sr0 ^ v {u-v} and U ref /* \J u ,veS r ,o,u^v{u- 

ct/N-: 
4V 2 



covers every nonzero element of GF(iV — 1) f (^-y- 2 - — 1) times. Thus the 
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q(N - l)/2 sets in {S r>a : r G U, a€ GF(N 
BIBD(/V - 1, q(N - l)/2, N/2 - 1). □ 



1)} constitute the blocks of a 



Proof of Theorem 2.2. Since T\ = ■■■ = T q / e and all the e sets in 
each Tj are distinct, where 7} = {T + (j — l)e, . . . , T + (j — l)e + e — 1}, 
1 < j ' < the e blocks {SVo : r = 0, . . . , e — 1} are themselves a difference 
family and are all distinct. Assume that certain two blocks of the design 
{S r ,a :r = 0, ...,e — 1, a G GF(iV — 1)} are the same. Then there exist an 
integer < r < e — 1 and a G GF(iV — 1), either r ^ or a ^ 0, such that the 
two sets Ki = {x™ +i + a : i G T, < j < (TV - 2)/g - 1} and K 2 = : i G 

T + r, < j < (iV — 2)/q — 1} are the same. Then the sum of the elements of 
K\ is equal to that of the elements of K2. Since q 7^ N — 2, 1 — x q 7^ 0; thus 
we have 

-(JV-2)/<Z-l 



(JV-2)/ ff - 

E 

J=0 



E 

j=0 



E^ 



0. 



This implies that the sum of the elements of K\ is equal to \{N — 2) a 
and the sum of the elements of K2 is equal to 0. Therefore, ^(N — 2) a = 
in GF(N — 1), and hence a = 0. Now since K\ = K2, the two sets K\ = 
{jq + i : i G T, < j < (TV - 2)/g - 1} and i? 2 = {j'q + i : i G T + r, < j < 
(N — 2)/q — 1}, with the elements being integers modulo N — 2, must be 
equal. By the definition of e, this can happen only if r = 0, which is a 
contradiction. □ 



Theorem 2.3 can be proved in the same way as Theorem 2.2. 

PROOF of Theorem 3.1. First we state and prove a lemma. 

Lemma 1. Let y be a 1 x (N — 1) vector with integer entries such that 
the first p— 1 entries are congruent to 2 (mod 4), and the last N —p entries 
are multiples of 4. Suppose p < N/2 and m' is the sum of the entries of 
y. If N - 1 < \m'\ < 2(N - 1) and 6 = (\m'\ -2p + 2)/4 is an integer, then 
< 9 < N — p, and the sum of squares of the entries ofy is minimized if and 
only ify has p — 1 entries equal to —2, 6 entries equal to —4 and N — p — 9 
entries equal to when m' < 0, or p — 1 entries equal to 2, 9 entries equal 
to 4 and N — p — 9 entries equal to when m! > 0. 

Proof. It is enough to prove the case m! < 0. The other case follows by 
reversing the signs. 
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If ml < 0, then 9 = (-m'-2p + 2)/4 < (2N - 2 - 2p + 2)/4 <N-p. Also, 
\m'\ > N - 1 6>[N-l-2p + 2]/4 > [N - 1 - N + 2]/4 > 0. Therefore, 
is a positive integer less than N — p. If a vector has p—1 entries equal to 
—2, 9 entries equal to —4 and N — p — 9 entries equal to 0, then the sum of 
all the entries is — A9 — 2(p — 1) = m! . We shall show that a vector y* with 
the smallest sum of squares of the entries among all vectors satisfying the 
conditions in the lemma must be of this form. 

First we show that y* cannot have positive entries. If not, let y* be a 
positive entry of y*. We claim that y* has at least one entry, say yj, such 
that i/j < —4. For otherwise, all the negative entries of y* are greater than 
or equal to —2. Then since p < N/2, at most N/2 entries of y* can be —2. It 
follows that ml - y* > (-2)N/2 = -N, but ml < -N + 1 and y* > 2 imply 
that ml — y * < —N — 1, a contradiction. Therefore, there is at least one 
y*A such that y* < —4. Now replacing y* and y* with y* — 4 and y* + 4, 
respectively, keeps the sum of the entries of the vector unchanged, but (y* — 
4 ) 2 + (Vj + 4 ) 2 < (V*) 2 + (Vj) 2 since \y* - 4| < \y*\ and \y* + 4| < \y*\. This 
means that y* can be improved, contradicting the fact that y*is optimal. 
Therefore, y* cannot have positive entries. 

Finally we show that y* can have only entries from the set {—4, —2, 0}. Let 
y*,..., y*N-i be the entries of y*. From the previous paragraph we know that 
all the entries of y* are nonpositive. Thus y* < —2 for all 1 < i < p — 1 and 
y j < for all p < j < N — 1. We first show that y* = —2 for all 1 < % < p — 1. 
Compare y* with the vector that has p — 1 entries equal to —2, 9 entries 
equal to —4 and N —p — 9 entries equal to 0. Since the sum of all the entries 
is a constant, we see that if there is a y* < —2, 1 < % <p — 1, then there must 
be at least one y* = 0, where p < j < N - 1. Then (y* + 4) 2 + (y* - 4) 2 = 
(y*) 2 + 8y* + 32 < (y*) 2 + (y*) 2 . The last inequality follows from y* < -6. 
This again shows that y* can be improved, which is not possible. Thus we 
must have y\ = ■ ■ ■ = yp_i = —2. Then the minimum of Y^i^iut) 2 subject 
to the constraint that all the y*'s are multiples of 4 is attained when each y* 
is or -4, p < i < N - 1. Since Y.fji y*i = m' , 9 = (-m' -2p + 2) /4 entries 
must be equal to —4. □ 

Now we are ready to prove the theorem. We denote the sum of squares 
of all the entries of a matrix M by SS(M). Then for a supersaturated de- 
sign X with m factors and N runs, E(s 2 ) = [SS(X r X) — mN 2 ]/[m(m — 1)]. 
A key fact used in Nguyen (1996) and Cheng (1997) is that SS(X T X) = 
tr[X r XX T X] = tr[XX r XX r ] = SS(XX r ), and since each column of X has 
the same number of l's and — l's, XX T has zero row sums. 

If N = 2 (mod 4), then all the entries of X r X are congruent to 2 (mod 4). 
In particular, all the off-diagonal entries have absolute values at least 2. 
Therefore, we have the simple lower bound 

(4.1) iV = 2(mod4) ^> £(s 2 )>4. 
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Let qi be the integer such that N — 1 < m — q\(N — 1) < 2N — 2. Since 
exactly one of m + q\, m + q\ + 1, m + q± + 2 and m + (/i + 3 is congruent to 
2 (mod 4), there is a unique q such that -2iV + 2 < ?n - q(N - 1) < 2iV - 2 
and m + q = 2 (mod 4). Let X* be obtained by adding q columns of l's to 
X and let J be the N x N matrix of l's. Then 

(4.2) X*(X*) T = XX T + gJ, 

and so SS(X*(X*) T ) = SS(XX T ) + g 2 SS(J) + 2q • (the sum of all entries of 
XX T ) = SS(XX T ) + q 2 N 2 . The last equality follows from the fact that 
XX T has zero row sums. Thus 

(4.3) E(s 2 ) = [SS(X*(X*) T ) - q 2 N 2 - mN 2 ]/[m(m - 1)]. 

A lower bound for E(s 2 ) can be obtained by bounding SS(X*(X*) T ). 

Without loss of generality, assume that each of the first p rows of X* has 
an even number of entries equal to 1 and each of the last N — p rows of X* 
has an odd number of entries equal to 1. We can also assume that p < N/2, 
since if needed we can change the signs of all the entries in a certain column 
of X*. Then X*(X*) T has the form 



(4.4) 



A C 
C T B 



where A is p x p, all the entries of C are multiples of 4 and all the entries of 
A and B are congruent to 2 (mod 4). This follows from the fact that m + q, 
the number of columns of X*(X*) T , is congruent to 2 (mod 4). 

Since XX r has zero row sums and its diagonal entries are equal to m, the 
sum of the off-diagonal entries in each of its rows is —m. Therefore, by (4.2) 
the sum of the off-diagonal entries in each row of X*(X*) T is —m + q(N — 1), 
which is a multiple of 4 if and only if q is odd and N = 2 (mod 4), since 
m + q = 2 (mod 4). Then since all the entries of C in (4.4) are multiples of 
4, 2(p — 1) is a multiple of 4 if and only if q is odd and N = 2 (mod 4). It 
follows that 

(4.5) p is odd if q is odd and N = 2 (mod 4); otherwise, p is even. 

Now we first consider the case where \m — q(N — 1)| < N — 1. Since all 
the entries of A and B are congruent to 2 (mod 4) , where A and B are as 
in (4.4), they all have absolute values at least 2. So SS(X*(X*) T ) is at least 
(m + q) 2 N + F(p), where 

F{p) = 4p(p - 1) + 4(7V - p)(N - p - 1) = 8p 2 - 8Np + AN 2 - AN. 

Since F'(p) = 16p — 8N has a zero at p = N/2 and F(p) is a convex function 
of p, by (4.5), F(p) is minimized at p = N/2 if N is a multiple of 4 or if N = 2 
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(mod 4) and q is odd, and at p = N/2 — 1 when N = 2 (mod 4) and q is even. 
From these observations we can calculate a lower bound for SS(X*(X*) r ): 



The various lower bounds for the case \m — q(N — 1)| < N — 1 as stated in 
the theorem can be obtained by combining (4.6) with (4.1) and (4.3). 

Next we consider the case \m — q(N — 1)| > N — 1. By the discussion in the 
paragraph preceding (4.5), the sum of the off-diagonal entries in each row of 
X*(X*) T is -m + q(N-l) and (\-m + q(N- 1)| - 2p + 2)/4 is an integer. By 
Lemma 4.1, with m! = —m + q(N — 1), the sum of squares of the off-diagonal 
entries of the first p rows of X*(X*) T is minimized if in each of these rows, 
p — 1 entries have absolute values equal to 2, (|— m + q(N — 1)| — 2p + 2)/4 
entries have absolute values equal to 4, and the rest are equal to 0. Thus 
SS(X*(X*) T ) is at least (m + q) 2 N + F(p), where 

F(p) = 4p(p- 1) + 4(N -p)(N -p - 1) 

+ 2 • (-4) 2 • p(\-m + q(N - 1)| - 2p + 2)/4 
= -8p 2 -8Np + 4N 2 - 4N + 8p(\-m + q(N - 1)|) + 16p. 

Since F'(p) = -8(2p + N — \—m + <?(iV - 1)| - 2) has a zero at p = [|-m + 
<7(iV — 1)| — N + 2]/2 and F is a concave function of p, by (4.5), F(p) is 
minimized at or N/2 [when iV = (mod 4)], or N/2 - 1 [when N = 2 
(mod 4) and q is even], and 1 or N/2 [when N = 2 (mod 4) and q is odd]. 
Here we have used the fact that < (|-m + g(7V - 1)| - N + 2)/2 < N/2. 

When TV = (mod 4), F(p) is minimized at if (|-m + q(N - 1)| - N + 
2)/2 > N/4, that is, if \—m + q(N - 1)| > |iV - 2; otherwise, it is minimized 
at p = N/2. Similarly, when N = 2 (mod 4) and q is even, F(p) is minimized 
at if |— m + q(N — 1)| > |iV — 3; otherwise, it is minimized at p = N/2 — 
1. When N = 2 (mod 4) and q is odd, F(p) is minimized at 1 if |— m + 
<7(iV — 1)| > |iV — 1; otherwise, it is minimized at p = N/2. Lower bounds 
for SS(X*(X*) T ) based on these observations together with (4.1) and (4.3) 
establish the various lower bounds for the case \m — q(N — 1)| > N — 1 as 
stated in the theorem. □ 
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